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We study Coulomb blockade oscillations of thermoelectric coefficients of a single electron tran- 
sistor based on a quantum dot strongly coupled to one of the leads by a quantum point contact. 
At temperatures below the charging energy Ec the transport of electrons is dominated by strong 
inelastic cotunneling. In this regime we find analytic expressions for the thermopower as a function 
of temperature T and the reflection amplitude r in the contact. In the case when the electron spins 
are polarized by a strong external magnetic field, the thermopower shows sinusoidal oscillations as 
a function of the gate voltage with the amplitude of the order of e~^|r|-^. We obtain qualitatively 
different results in the absence of the magnetic field. At temperatures between Ec and _Ec'|''P the 



thermopower oscillations are sinusoidal with the amplitude of order e 



' In . On the other 



hand, at T ^ Ec\r\ we find non-sinusoidal oscillations of the thermopower with the amplitude 

^e'^\ryT/Ec\n{Ec/T). 

PACS numbers: 73.23.Hk, 73.50.Lw, 72.15.Jf 



I. INTRODUCTION 

It is well known that electric current in solids can be 
caused not only by an applied electric field, but also by 
a temperature gradient. This gives risft,to a number of 
interesting thermoelectric phenomena.BH It is important 
to note that the thermoelectric phenomena in metals re- 
quire asymmetry between electrons and holes. Indeed, 
in a perfectly electron-hole symmetric system the tem- 
perature gradient will cause equal in magnitude currents 
of electrons and holes which result in zero net electric 
current. 

Many recent studies of thermpelectric effects were fo- 
cused on mesoscopic systems. utJ Thermoelectric prop- 
erties of these systems are particularly interesting be- 
cause the electron-hole asymmetry in mesoscopic devices 
is usually strong and can be controlled experimentally 
by tuning external parameters, such as gate voltage or 
magnetic field. 

Of the various thermoelectric phenomena the Peltier 
effect is probably the most important for technological 
applications: when an electric current / is passed through 
a system in the absence of the temperature gradient, it 
is accompanied by the heat current 
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Here 11 is the Peltier coefficient. The use of the Peltier 
effect has been proposed for refrigeration in conditions 
when various technological constraints, such as the size of 
the device, outweigh the power efficiency considerations. 
The strong enhancement of the particle-hole asymmetry 
in mesoscopic devices and their small size make them 
very promising candidates for microrefrigerators.l3 



In the last Jew years many experimental and theoret- 
ical studieaTEj focused on the thermoelectric properties 
of single electron transistors (SET). Thermoelectric ef- 
fects in these systems can be controlled by the gate volt- 
age Vg, Fig. 1^, and exhibit characteristic Coulomb block- 
ade oscillations. Most of the studies of thermoelectric 
effects in the Coulomb blockade regime concentrated on 
the thermopower 
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Here V is the voltage induced across the device in the 
absence of net electric current when the temperatures of 
the two leads differ by AT, Fig. ^(a). The Peltier coef- 
ficient H is related to the thermopower S by an Onsager 
relation H = ST. 

The theory of the Coulomb blockade oscillations in the 
thermopower of single electron transistors in the weak 
tunneling regime was developed in Ref. This theory 
takes into account only the lowest order tunneling pro- 
cesses, i.e. the sequential tunneling, and neglects the co- 
tunneling processes. Its results ysere in agreement with 
the experiments of Ref. ^. LaterO it became possible to 
experimentally access the regime of lower temperatures 
and stronger coupling to the leads, where the cotunneling 
processes become dominant. The theoretical description 
of this regime was recently given in Ref. 

In experiments with GaAs heterostructures the quan- 
tum dot is connected to the leads by quantum point con- 
tacts. Each contact is usually in the regime when only 
one transverse mode can propagate through it, and the 
transmission coefficient fpi: this mode can be controlled 
experimentally. Recentlyll3 the Coulomb blockade oscil- 
lations in the thermopower of a SET with the quantum 
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dot strongly coupled to one of the leads were studied 
for various values of the reflection coefficient |rp in this 
contact. The setup of these experiments is schemati- 
cally shown in Fig. |l](b). In the regime of strong cou- 
pling, |rp ^ 1, nearly sinusoidal oscillations of the ther- 
mopower as a function of the gate voltage Vg were ob- 
served. 
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FIG. 1. (a) Setup of the thermopower measurement in a 
single electron transistor. A quantum dot is capacitively cou- 
pled to the gate and connected to the two leads by tunneling 
junctions. The left and right leads are maintained at temper- 
atures T + AT and T, respectively, and the voltage V across 
the device is measured. The thermopower (j^ is measured as 
a function of the gate voltage Vg. (b) A SET with a quan- 
tum dot strongly coupled to the right lead by a single-channel 
quantum point contact. The cross x in the constriction rep- 
resents the backscattering in the contact resulting in a finite 
reflection coefficient |rp. 

The theories of thermopower for the weak tunneling 
regime developed in Refs. 0,Qdlrely upon the perturbation 



theory in the strength of coupling between the quantum 
dot and the leads. This perturbative approach fails when 
the coupling to the leads is strong. The previous theo- 
retical work on Coulomb [blockade systems with strongly 
coupled quantum dotal3ila was devoted to the studies of 
their thermodynamic properties and conductance. These 
properties are not sensitive to the electron-hole asymme- 
try, and the calculation of the thermoelectric properties 
requires a non-trivial generalization of the approach of 
Refs. 

In this paper we develop a theory of the thermopower 
in a SET with a quantum dot strongly coupled to one of 
the leads, Fig. |l](b). We consider a relatively large quan- 
tum dot in which the quantum level spacing is small, and 
the electron transport is dominated by inelastic cotun- 
neling. In Sec. ^ we present a qualitative discussion of 
thermoelectric transport in a SET. The thermopower of a 
SET in the regime of strong inelastic cotunneling is devel- 
oped in Sec. Ill for the simpler case of spin-polarized elec- 
trons as well as for the more interesting spin-degenerate 
case. We discuss the results and compare them with ex- 
periments in Sec. IV. 



II. QUALITATIVE DISCUSSION OF THE 
THERMOPOWER AT WEAK INELASTIC 
COTUNNELING 

The physical meaning of the thermopower can be de- 
duced from the Onsager relation S = II/T. One can 
easily see from Eq. (|l]) that the Peltier coefficient is deter- 
mined by the average energy (e) of the electrons carrying 
current through the system, measured from the chemical 
potential: 11 — — (e)/e. Thus the thermopower measures 
the average energy of the tunneling electrons in units of 
the temperature: 



(3) 



Here e is the absolute value of the electron charge. The 
average energy of the charge carriers (e) is determined by 
a particular mechanism of transport through the system. 

A conventional SET schematically shown in Fig. |^(a) 
consists of a quantum dot weakly coupled to two leads. 
The transport of electrons from the left lead to the right 
one is achieved by either sequential tunneling or cotun- 
neling (elastic or inelastic). 

The sequential tunneling refers to the lowest-order tun- 
neling processes in which one electron tunnels into or out 
of the dot. As a result of each tunneling event the charge 
of the dot changes by ±e. When an electron tunnels into 
or out of the dot, the electrostatic energy of the system 
increases by u+ or u_ , respectively. The values of u+ and 
U- are of the order of the charging energy Ec = e^/2C, 
where C is the capacitance of the dot; their values can 
be tuned by adjusting the gate voltage Vg. The electron 
that tunnels into the dot has to have the energy e « in 
order to charge the dot, so the tunneling is exponentially 
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suppressed at low temperatures as e~"+/"^. Similarly the 
rate of tunneling out of the dot is suppressed as e~^^^^ . 
The thermopower of a SET in the regime of sequential 
tunneling was studied in Ref. 0. 

The cotunneling mechanism accounts for the coherent 
second-order tunneling processes in which at the first step 
an electron tunnels from, say, the left lead into the dot, 
and at the second step this or another electron (for elas- 
tic and inelastic cotunneling, respectively) tunnels from 
the dot to the right lead. Since as a result of cotunnel- 
ing processes the charge of the dot remains unchanged, 
the electrons participating in the transport through the 
dot do not need to have large energies ^ u±, and the 
transport is not suppressed exponentially at low temper- 
atures. It is important to note that the elastic cotunnel- 
ing involves elastic propagation of electrons between the 
tunneling contacts in the dot. The resulting contribution 
to the transport is inversely proportional to the volume 
of the dot. Thus in relatively large dots the finite tem- 
perature transport is dominated by inelastic cotunneling 
processes. 

In order to find the cotunneling thermopower of a SET, 
one needs to evaluate the average energy (e) of the elec- 
tron transferred through the device. Inelastic cotunnel- 
ing involves electrons within the energy strip of width 
~ T near the Fermi level. The probability w(t) of the 
second-order tunneling process is inversely proportional 
to the square of the energy of the virtual state. 
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where we assumed that the electron with energy e tun- 
nels into a state of energy e' in the dot. The energies e 
and e' are of the order of T and small compared to the 
charging energies u± . One can therefore expand Eq. (Q) 
in small (e — e')/u±, 
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The e-dependent correction in Eq. (||) shows that the 
tunneling probability increases or decreases with the en- 
ergy e of the tunneling particle depending on the sign of 
{u^^ — uZ^). Since the typical energy e ^ T, the relative 
magnitude of the term breaking the electron-hole sym- 



metry in (|^) is ~ T{i 



~_^) and the average energy 



of the tunneling particles is (e) '-^ T 



:^). Then 



using Eq. (H) we estimate the thermopower of a SET as 
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Here A is a numerical coefficient of order unity; its value 
A — 47r^/5 has been found in Ref. ^ 

An interesting feature of the result (||) is that the ther- 
mopower does not depend on the strength of coupling of 
the dot to the leads. This feature is expected to persist 



as long as the transmission coefficients of the barriers are 
small. However, as the barrier approaches the regime of 
perfect transmission, the Coulomb blockade osciLla.tions 
of physical quantities are expected to disappear .E3o One 
should therefore expect that in the regime of strong cou- 
pling the thermopower will depend on the reflection co- 
efficient |rp of the contact. 

It is also worth mentioning that the thermopower 
S ^ e^^T/Ec in a SET is much greater than the typ- 
ical value S ~ e~^T/ Ep of the thermopower in metals; 
here Ep \s the Fermi energy. This is a consequence of 
the fact that the charging effects in the dot enhance the 
electron-hole asymmetry. We will show in Sec. [II that 
such behavior of the thermopower at T — s- persists in 
the strong coupling regime. 



III. THERMOPOWER IN THE REGIME OF 
STRONG INELASTIC COTUNNELING 

In this section we calculate the thermopower S" of a 
SET in which one of the contacts is in the strong tun- 
neling regime. Fig. |^(b). We restrict ourselves to the 
linear response regime when both the temperature dif- 
ference between the leads AT and the voltage V across 
the device are small. The thermopower (||) is defined in 
terms of the voltage V induced across the device by the 
temperature difference AT at zero current. In practice, 
however, it is easier to calculate the current response 



I = GV + GtAT. 
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Here G is the conductance of the SET, and Gt is the 
thermoelectric coefficient describing the current response 
to an applied temperature difference. The thermopower 
(0) can then be expressed as 



Gt 
~G' 



(8) 



The conductance G of a SET in the regime of strong in- 
elastic cotunneling was found in Ref. W% Hence in the 
following we concentrate on the calculation of the ther- 
moelectric coefficient Gt- 



A. Tunneling approximation 

We assume that the conductance of the tunneling 
junction connecting the dot to the left lead is much 
smaller than the conductance quantum, G^ ^ /h. In 
this case one can describe the SET by the Hamiltonian 
H = Hl + H' , where Hl is the Hamiltonian describing 
the tunneling of electrons in the left contact. 
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and H' accounts for the transport through the right junc- 
tion and the electron-electron interactions in the system. 
Here au and Op are the electron annihilation operators 
in the left lead and in the dot, respectively, tk and 
are the energies of the corresponding states, the matrix 
elements t^p describe the weak tunneling of the electrons 
through the barrier. 

We will account for the tunneling through the left bar- 
rier in the lowest (second) order of the perturbation the- 
ory in tkp- In this approximation all of the temperature 
drop occurs at the tunneling barrier. We take the tem- 
perature of the left lead to be T-l- AT and that of the dot 
to be T. Let us denote the tunneling density of states in 
the left lead by vi{e) and that in the dot by v{e). To the 
lowest order in the tunneling matrix element the current 
through the tunneling contact can be obtained with the 
aid of the Fermi golden rule: 
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ui(e)v(e)[ni[e)-n{e)]de. (10) 



Here ni{e) and n(e) denote the Fermi distribution func- 
tions at the temperature of the left lead T + AT and the 
dot T, respectively. The square of the tunneling matrix 
element in ] 
Fermi level. 

To determine Gt we assume that the chemical poten- 
tials in the Fermi functions in Eq. ( p^ ) are the same and 
expand the current to first order in AT, 



veiy. 

element in Eq. (|l^) is averaged over the states near the 
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Here we have replaced the density of states in the left 
contact vi{e) by its value vi at the Fermi energy using 
its weak energy dependence. The corrections to this ap- 
proximation are small in the ratio of the temperature to 
the Fermi energy. We will see below that the density 
of states in the dot has the energy dependence at the 
much smaller energy scale T due to the electron-electron 
interactions. 

In order to express the thermoelectric coefficient Gt 
in terms of physically measurable quantities, we use 
Eq. ( p^ ) to calculate the conductance Gl of the left 
barrier assuming that the electrons in the dot are non- 
interacting. This conductance can, in principle, be mea- 
sured experimentally by opening completely the constric- 
tion connecting the dot to the right lead. The result is 
expressed in terms of the density of states in the dot vq, 
which is no longer renormalized by the electron-electron 
interactions: 



GL = 2'Ke'uiVa{\tkp\'')- 



(12) 



Using Eqs. (|T^) and ( |ll| ) we can now express Gt — I j 
as 
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Equation (|1^) reduces our problem to the calculation 
of the energy-dependent tunneling density of states J^(e). 
We note that unlike the conductance G of the SET, the 
thermoelectric coefficient Gt is determined by the odd 
(in energy) component of density of states i^(e). There- 
fore the thermopower measurements represent an inde- 
pendent test of the theory of Coulomb blockade in nearly 
open dots developed in Refs. 13-1^. 

It is well known that the tunneling density of states can 
be expressed in terms of the electron Green's function. A 
specific form of this relation that will be convenient for 
further calculations is 
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Here g{T) = -(TrV'L(T)V'l(O)) is the Matsubara Green's 
function; tpL the annihilation operator of an electron in 
the dot at the position of the left contact. For the deriva- 
tion of E n. (H ) see Appendix |^. Substituting Eq. (|l|) 
into Eq. (13|), we express the thermoelectric coefficient 
Gt in terms of the Green's function: 

Gt = / T G { — +it]dt. (15 

2et/o 7-00 cosh^ {nTt) \2T J ^ ^ 

This expression is insensitive to the specific form of the 
interactions in the quantum dot and its coupling to the 
right lead. In the following sections we calculate the 
Green's function C/(t) in the strong inelastic cotunnel- 
ing approximation and find the corresponding value of 
Gt- 



B. Inelastic cotunneling approximation 

At finite temperature in a sufficiently large dot one can 
use the approximation of inelastic cotunneling which ne- 
glects the possibility of elastic propagation of electrons 
between the two contacts in the dot. In this case the 
system can be modeled by the Hamiltonian 



H — Hl + Hr + He , 



(16) 



in which the Hamiltonians and Hr describe the in- 
dependent subsystems of electrons propagating through 
the left and right contacts; in particular [Hl , Hr] = 0. 
We have discussed the form of H^ in section III A, see 
Eq. For weak inelastic cotunneling Hr would have 
a form similar to Eq. (||). In the case of strong inelastic 
cotunneling the Hamiltonian Hr has a comp letely dif- 



ferent form which will be discussed in sections III C and 
Finally, the term He describes the Coulomb inter- 
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actions in the dot. At low energies the interactions are 
adequately accounted for by the charging energy approx- 
imation: 



He = EcinL + fiR- 



(17) 
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Here fi^ and fiR are the operators of the number of elec- 
trons that entered the dot through the left and right con- 
tacts, respectively, and is a dimensionless parameter 
which is proportional to the gate voltage Vg. 

By definition of the operator its commutation rela- 
tions with the fermion operator ij^L defined below Eq. jl^ ) 
have the form [-0/, ^fi^] = ipL ■ For the convenience of the 
following calculations we will rewrite the charging energy 
fO) in the form 



Hc = Ec{n + nR-N)^ 



(18) 



where n is an integer-valued operator that commutes 
with ipL. In order to preserve the commutation rela- 
tions between V'l and He we replace t/'l '4'lF, where 
F is the operator lowering n by unity: [F, n] = F. It 
is important to note that the substitution t/'l ^lF 
does not affect the form of the Hamiltonia n Hl , and the 
only modification in the discussion of Sec . [Ill A | is in the 
definition of the Green's function in Eq. (^4| 
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The operators ipL and now commute with Hn + He, 
whereas F and i^^ commute with H^. Consequently, the 
Green's function ( p^ ) factorizes, 

Gir) = -(r.VL(r)Vl(0))(r.F(r)Ft(0)). 

In the new representation the operators ipL and de- 
scribe non-interacting fermions, whose Green's function 
is well known. We can then rewrite G{t) as 



vqttT 
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sin(7rrT) 
if(r) = (T.F(r)Ft(0)). 



(20) 
(21) 



Substituting Eq. (|20|) into Eq. (|15|) we express the ther- 
moelectric coefficient Gt of the dot in terms of the cor- 
relator K{t), 



Gt = — 
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sinh TrTi) / 1 , , 
'-K \-it]dt. 22 



Unlike Eq. (|T|), formula (||) assumes that the transport 
through the dot is due to the inelastic cotunneling mech- 
anism, and that the electron-electron interactions in the 
system are completely described by the charging energy 
(p^). On the other hand, the coupling Hj^ of the dot to 
the right lead is still arbitrary. In particular, Eq. (|2^ ) is 
valid in the case of a metallic grain coupled to the lead 
by a wide contact supporting many channels. 

In this paper we consider the case of a single-channel 
contact, which is usually realized in semiconductor de- 
vices. Depending on the presence of a magnetic field 
polarizing the spins of the electrons, one has to consider 
the cases of either spinless or spin-i electrons. 



C. Spinless electrons 

We start with the simpler case of spinless electrons. 
Following Ref. |l^,|l^ we describe the electron transport 
through the right quantum point contact by a model of 
one-dimensional fermions. In the case of strong cou- 
pling of the dot to the right lead, the charging energy 
( |l8| ) gives rise to non-trivial Coulomb correlations of 
the motion of electrons through the constriction. It is 
more convenient to treat the problems of interacting one- 
dimensional electrons-iii—the bosonized representation. 
Then the HamiltoniarOa of the right constriction takes 
the form Hr = ijjj^-* + i/^, where 

/oo 
{TT^n^{x) + [d,cf>{x)]^}dx, (23a) 
-oo 
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R 



r\cos[2(t>{0)]. 



(23b) 



Here (j) and H are bosonic fields satisfying the commu- 
tation relations [(/)(a;), n(?/)] = iS{x — y), the parameter 
Ui? is the Fermi velocity of the electrons, r is the reflec- 
tion amplitude in the constriction, and D is the band- 
width. The regions x < and a; > in the integral of 



Eq. (23a) represent the electrons in the dot and in the 
right lead, respectively. The deviation of the density of 
one-dimensional electrons from its ground state value is 
given by dx4>{x)/T:. Thus the number of electrons that 
have entered the dot through the right constriction is 
"fl — 0(0)/^: and the charging energy ( psf ) takes the 
form 



Hc = Ec 



1 



(23c) 



The advantage of the bosonization approach is that 
the Coulomb interaction term (23c) is quadratic in the 
bosonic operator 0, and, therefore, can be treated ex- 
actly. On the other hand, the backscattering of elec- 
trons in the constriction in the bosonized representation 
takes the strongly non- linear form ( ^3b| ) . As a result the 
backscattering can only be accounted for perturbatively, 
using the small parameter \r\ < 1. 

We will calculate the time-ordered correlator K{t) de- 
fined by Eq. (j^) as an imaginary-time functional integral 
over the bosonic field <j). The operator i^^(O) increases n 
from to 1 at time < = 0, whereas F{t) changes it back 
to n = at time t = t. Therefore i^(0)F'l'(0) in the func- 
tional integral can be omitted provided that the operator 
n in the action is replaced by 



ririt) ^9{t)0{T-t). 



(24) 



Here 9(t) is the unit step function. Upon this procedure 
the correlator K{t) is expressed as follows 



Zjr) 

z{0)- 



(25) 
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where Z{t) is the functional integral given by 



Z(r) 



exp[-5o - Sc{t) - S']V(j>{x, t). 



(26) 



Here Sq denotes the part of th e Eucledian action de- 
rived from the Hamiltonian (23a) of free electrons moving 
through the constriction in the absence of both interac- 
tions and backscattering: 



So 



dt I dx— ^ 
27r 



(27a) 



Here (3 = 1/T. The term Sc{t) is the par t of the ac- 
tion which is due to the charging energy (23c), where the 
operator h is replaced with (P4), 



Sc{r) 



Er 



nr{t) + -(t>{Q,t)-N 

TT 



dt. (27b) 



Finally, S' is the s mall c ontribution to the action due to 
the backscattering (23b), 



S' 



D 



r|cos[2(/>(0,t)]dt. 



(27c) 



The following calculations are performed in the regime 
of low temperatures, T <C Ec- At r = the functional 
integral (E6) is gaussian, and its explicit evaluation gives 



^o(t) 



1 



\lEc J sin2(7rrr)' 



(28) 



see Appendix yj. Here 7 = e*^, where C « 0.577 is 
Euler's constant. The substitution of ( p8| ) into the ex- 
pression ( p2| ) results in Gt = 0, as the integrand is an 
odd function of t. This is a consequence of the fact 
that at r = Othe system described by the Hamiltonian 
( p3a| ) and (23c) possesses electron- hole symmetry. In- 
deed, in the absence of backscattering the dependence of 
the Hamiltonian on the gate voltage N can be removed 
by shifting the field <t){x) (j){x) +ttN. Then the Hamil- 



tonian H 



He is obviously invariant with respect 



to the electron-hole symmetry transformation n — > —h, 
(j) — > —(f), and H — H. On the other hand, the shift 
-> 0(x) + ttN will change the backscattering term 



( 231: ) . Thus the backscattering breaks the electron-hole 
symmetry and gives rise to non-zero thermopower of the 
device. 

To account for the small backscattering at |r| <C 1 we 
expand the action in Eqs. (^5|) and (|2^) to first order in 
S' and find 



K{T)=Ko{T){l-{S')r + {S')o). 

Here {S')r is defined as 



{S')r 



J exp[-5o-5c(r)]I?0 



(29) 



(30) 



Due to the form (27c) of the perturbation S' the eval- 
uation of K{t) in Eq. ( p9| ) again amounts to taking 
gaussian functional integrals. The straightforward but 
lengthy calculations carried out in Appendix ^ give the 
result 



Kir) 



Xo(r)[l - 27^|r| cos(27riV) 

-f47r2^7|r| ^ sin(27riV) cot(7rrr)] . (31) 



Here ^ « 1.59 is a constant defined by Eq. ( p31l| ). The 
substitution of this result into Eq. (|2^) gives the following 
result for the thermoelectric coefficient in the first order 



(32) 



It is interesting to note that the thermoelectric coeffi- 
cient Gt given by Eq. (|3^) is an odd function of the gate 
voltage N. This property is more general than the per- 
turbative result ( |3^ ) . Indeed, one can see from the form 
of the functional integral ( |26| ) and the action ( p7| ) that 
the correlator (E3) has the following symmetry property: 
K{P-T, N) = k(r, 1-iV). Furthermore, aU the physical 
properties of the system are periodic in N with period 1. 
This can be shown by shifting — s- + nN which re- 
moves N from Sc (t) and changes the cosine in Eq. ( p7c| ) 
to cos[20(O, t) + 2TTN]. The action then becomes invariant 
with respect to the shift N ^ N + 1. Consequently the 
correlator ( |25| ) has the property K{P—t, N) — K{t, —N). 
One can easily see from Eq. (22) that only the part of 
K{t) which is odd with respect to r ^ P — t contributes 
to the thermoelectric coefficient Gt- Using the afore- 
mentioned properties of K{t,N), this odd part can be 
presented as 



Koddir) - (r, N) - K{p - r, N)] 



[K{t,N)-K{t,-N)]. 



(33) 



Therefore the thermoelectric coefficient Gt is an odd 
function of N. 

The conductance G can be obtained by substitution 
of correlator ( ^l|) an alytically continued to real time into 
Eq. (52) of Ref. Wl The result has the form 



27r'*r2 



27^|r|cos(27riV)]. 



(34) 



This expression is in agreement with the formula (A27) 
of Ref. 1^, where the numerical prefactor in the brack- 
ets was not determined, and with the expression (34) of 
Ref. |l^, where the constant ^ was found. 

Substituting En. (^) into Eq. (^ and using the lead- 
ing term in Eq. ( |34| ) for the conductance, we obtain the 
following expression for the thermopower in the spinless 
case 



6 



S ■ 



' ^eEc 



r|sin(27riV). 



(35) 



It is instructive to compare this result with the ther- 
mopower (^ in the regime of weak inelastic cotunneling. 
Both expressions vanish linearly at T — > 0, but unlike 
Eq. (^), our result ( ^ ) depends on the transmission co- 
efficient of the barrier. As expected, at perfect trans- 
mission r ^ the Coulomb blockade oscillations of the 
thermopower disappear. 



D. Electrons with spin 

Although the spins of electrons can be polarized in an 
experiment by applying a strong magnetic field, the most 
common situation is when no field is applied. In this 
regime one has to consider the case of spin-i electrons. 

In the presence of electron spins the Hamiltonian 
( ^3a| )-( p3^ ) has to be modified to account for the two 
species of electrons: spin-t and spin-|. Each of the 
spin subsystems can be bosaoized independently, and the 
Hamiltonian takes the forrr£3 



H 



(0) 

' 2tt 



I {7^^K{^) + [dMx)f}dx, (36a) 

. I J —OQ 



H'r = \r\ {cos [20t(O)] + cos [2^0)]} , (36b) 

vr 



(36c) 



To find the thermoelectric coefficient (^) one has to find 
the correlator K(t). Similarly to the case of spinless elec- 
trons, K{t) can be expressed in terms of the imaginary- 
time functional integral (p5[), where 



Z{t) = j exp[-5o - Sc{t) - S']Vcj),Vcbs 



(37) 



Here we have introduced the charge and spin fields 
4>c,s{x^t) = [0|(x,t) ± 01 (x, i)]/V2. The action in 
Eq. ( ^ ) is expressed in terms of these variables as 



So = / dt dx 



VF 

2tt 



{dt(t>a 



(38a) 



Scir) 
5' = - 



dtEc 







/2 

TT 



H 2 



at \r \ cos 

I" 



(38b) 

\/20,(O,<)]cos[V20,(O,i)]. (38c) 



Similar to the case of spinless electrons, in the absence 
of backscattering in the constriction the calculation of 
K(t) reduces to evaluation of a gaussian functional in- 
tegral (|^. Clearly, at r = the integral over 0^ is 
unaffected by nr{t); therefore the integrals over the spin 
degrees of freedom in the numerator and denominator 



of E q. (25) can cel. One can easily see that the action 
( |38a| ) and (38b) of the charge mo de is identical to that 
of the spinless problem ( ^7aD and ( ^7b| ) upon the substi- 
tution Ec 2Ec, nr{t) nr{t)/V2, and N N/V2. 
Making the respecti ve m odifications to the derivation of 
Kq(t) in Appendix B 1, we find 



Kc(r) 



TT^T 



2fEc |sin(7rTT)| 



(39) 



Substituting the analytic continuation of this result to 
T = 1/2T + it into Eq. ( ^, we find Gt = 0. As it 



was explained in Sec. [II C, this is the consequence of the 
fact that the system possesses electron-hole symmetry at 
r = 0. 

Th e rest of this section is organized as follows. In 
Sec. [II D 1 we calculate the thermopower within the 
second-order perturbation theory in the reflection am- 
plitude r. We show that the perturbative result diverges 
at low temperatures. We then find the the rmopower at 
arbitrarily low temperatures in Sec. HID 2 using a non- 
perturbative approach. 



1. Perturbation theory 



At non-vanishing backscattering the correction to 
Kc{t) appears in the second order in r. Indeed, the first- 
order correction can be expressed in the form (p^). It is 
easy to check that unlike the case of spinless electrons, 
the average {S') vanishes, because the fiuctuations of the 
spin mode 0^(0, i) are not suppressed at low frequencies 
by the charging energy term (38b). Expanding Eqs. (|2j 
and (^) to second order in r, we find 



Kir) = K.ir) 



(40) 



Similarly to Eq. (|30|), the averaging {■■■)t he re is per- 
formed with the action Sq + S c (t) given by (38a) and 
( |38b| ). Using the explicit form (|38c[ ) of S', we get 



4^2|^|2 rP rl3 



Kc{t,t';T)Ks{t,t')dtdt', (41) 



'0 Jo 

where we have introduced the correlators 

Kjt, t'; t) = (cos[%/20c(O, t)] cos[%/20c(O, i')]>r , (42a) 
Ks{t,t') = {cos[V2(f>s{0,t)]cos[V2(f>s{0,t')])o. (42b) 

The spin fiuctuations are completely decoupled from the 
charging action (38b), rendering the correlator Ks inde- 
pendent of T. The calculation of the correlators Kc and 
Ks reduces to evaluation of gaussian integrals. In Ap- 
pendix C 1 we find 
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ttD 



K,{t,t') 



nT 



1 



2D\sm[TTT{t-t')]\' 
where we have introduced the notation 

Xrit) = TTUrit) + SXrit), 

sin[27rT(t — T)n] — sui[2T:Ttn\ 



Ec 



(43) 
(44) 

(45a) 
(45b) 



we have assumed |t — t'| S> ^^p^; we wih 



In Eq. 

see that this region gives the leading contribution to the 
integral (^). 



As discussed in Sec. [II C, only the odd in N part 



of the correlator K{t) contributes to the thermoelectric 
coefficient (p2[). Keeping only the odd part of Eq. (43), 
from (^) in the second order in r we find 

i^odd(T) = X,(r)^l^|r|2 sin(27r7V)X(T), (46) 

Jo Jo \sm[-KT{t-t')]\ 

In evaluating the integral T(r) one should keep in mind 
that the denominator in Eq. ( |47| ) is written for |i — 1'| ^ 
E^^. Thus the logarithmic divergence at i = t' should 
be cut off at \t-t'\ ~ 

To evaluate the integral X(t) we first notice that away 
from the points t = 0,t, /? the correction Sxrit) in 
Eq. (H) is small in T/Ec, 



6Xr{t) 



2E^ 



{cot[7rr(t - r)] - cot[7rrt]}. 



(48) 



One can neglect this correction in the argument of the 
cosine in Eq. (^ and replace cosxr(t') = sgn(i' — r). 
Then the integral over t' can be evaluated with logarith- 
mic accuracy: 



2 ff^ 
- J d,t sin Xrit) 



— In tan 



sgn(< - t) In ( ^ tan ^^^^^ '^^ 



To leading order in T/Ec one can replace sinxT(i) — 
Sxrit) sgn(t — r). Using the approximation ( ]4^ ) we then 
obtain with logarithmic accuracy 



i^odd(r) = -A'c(r)§|r|2sin(27rAr) In In tan 



ttTt 



(49) 



Substituting this result for K{t) in Eq. ( p2| ) we find 
9 e ^ T 



Gj, = _^^^ln:^|rpsin(27r7V). (50) 



The second-order perturbation theory result for the ther- 
mopower can be found from Eq. (|^) using the result 

^ 8fEc 

1 



G^G 
at r 



of Ref. 14 for the conductance of the device 



^--S-ln^M^sin(2.iV). 



97r2 



T 



(51) 



This result applies at T ^ Ec and, similarly to the spin- 
less case (^5|), the thermopower vanishes at r ^ 0. It is 
important to note that unlike the spinless case (|35|), the 
thermopower ( |5l| ) diverges at T ^ 0. This means, in par- 
ticular, that the perturbation theory leading to Eq. (|5l| ) 
fails at sufficiently low temperatures. In the next section 
we perform a non-perturbative calculation and establish 
the true behavior of the thermopower at T — > 0. 



2. Non-perturbative treatment 



The logarithmic growth of the thermopower (^ij) at 
low temperature indicates that the thermoelectric prop- 
erties of the system are controlled by the spin and charge 
fluctuations at frequencies below Ec- In this section we 
construct a theory that describes the low-energy prop- 
erties of the system exactly and enables us to obtain a 
non-perturbative expression for the thermopower at ar- 
bitrarily low temperatures. This derivation was outlined 
in Ref. ^ 

As we already discussed, at r = the contributions 
of the spin fluctuations to the functional integrals in 
the numerator and the denominator of Eq. (|2^) cancel 
each other, and the ratio of the functional integrals over 
the charge degrees of freedom is equal to the correlator 
Kc{t), Eq. (^. The effect of smaU but finite r on the 
charge modes is negligible, because their fiuctuations at 
low energies are suppressed by the charging energy. How- 
ever, even a small backscattering r pins the fluctuations 
of the spin modes apd changes their low-frequency dy- 
namics dramaticallylia Therefore one can account for the 
small backscattering by presenting the correlator ( p5|) in 
the form 



K{t)^K,{t)Ks{t), K,{t) 



Zs{t) 



(52) 



where Zgij) is the functional integral over the slow spin 
modes, averaged over the fast charge modes. 

The calculation of Zs(t) amounts to integrating out 
the fast charge degrees of freedom in the functional in- 
tegral (|37|). Since the spin and charge fl uctu ations are 
only coupled by the backscattering term (38c), this pro- 
cedure reduces to the averaging of cos[-\/20c(O, t)] with 
the gaussian action So + Sc(t). Indeed, one can rewrite 
Eqs. dH) and (B|) as 
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(cos[v^0,(O, t)] cos[y2<^,(0, t')])r = Tr{t)Tr{t'), (53) 

(54) 



One can se e fr om Eq. (^) that at |t — i'| ^ -E^^ the 
correlator (42a) factorizes into the product of the aver- 
ages of the cosines, and that (cos[-\/20c(O, t)])T = '^rit)- 
It is clear that the higher-point correlators will also fac- 
torize into the product of averages. One can therefore 
simply replace the cos[\/20c(O, i)] in the action ( ^8c| ) by 
Tj- (t) and obtain the effective action for the spin degrees 
of freedom in the form 




Xr{t)cos[V2(t)s{0,t)]dt, (55) 



/O V VF 

where we have introduced the notation 

Xrit) = AcOslXrit) - TtN] 

= A(-l)"-(*) coa[Sxrit) ~ ttN], 



A = 



\r\. 



(56) 
(57) 



The procedure leading to the action (^5|) implied that all 
the relevant time scales of the problem are longer than 
E^'^. Therefore one has to integrate out the fluctuations 
of the spin degrees of freedom with frequencies exceed- 
ing Ec- This procedure is straightforward and amounts 
to replacing D with the new bandwidth ~ Ec- Thus 
Eq. (^^ gives the effective action of the problem, pro- 
vided the bandwidth D ~ Ec- 

One can now find the correlator K{t) using Eq. ( |5^ ) 
where the functional integral Zs{t) is defined as Zs(t) = 
/ e~^^'D(l)s- For the subsequent calculations it will be 
convenient to use the hamiltonian formulation of the 
problem and express Zs{t) as the trace of the time- 
ordered exponential: 





r0 




Zs(t) = Tr |rt exp 


- / Hrit)dt 


1 , (58) 




Jo 





where the time dependent Hamiltonian is given by 

/oo 
-oo 




\r{t) C0s[V2(j)s{Q)]- 



(59) 



The small parameter of the problem r enters through 
Xr{t)- In order to evaluate K{t) in all orders in Ar we 
refermionize the Hamiltonian (^) following Ref. ^ and 
find 



Hrit) 



^kclck - Xr{tKc + c^){ck - cl)] dk- (60) 



Here = vpk; the operators c|, and Ck satisfying the an- 
ticommutation relations {cfc,cj.,} = 6{k — k') create and 
destroy chiral fermions. Finally, c is a fermion annihila- 
tion operator anticommuting with Ck and c^. 

Although the Hamiltonian ( |60|) is quadratic in the 
fermion operators, the time dependence of Xr(t) makes 
the evaluation of the trace (pSh non-trivial. It is clear 



from Eq. (45b) that at T/ Ec ^ 1 the main time depen- 
dence is due to the factor (— 1)"t(*) in the definition of 
Ar(t), Eq. (IsgI). One can greatly simplify the calculation 
by eliminating this time dependence with the following 
trick. Note that the unitary transformation with the op- 
erator 



U 



f)(c- 



ct) 



(61) 



changes the sign of At in the Hamiltonian (^0|). There- 
fore the factor (— l)"T(t) g^n be accounted for by adding 
operators U{t) and C/(0) to the trace (pi). 



Zs{t) = Tr < Tt exp 



[Hq ^ H'^{t)]dt 



U{t)U{Q) 



(62) 

Here HQ+H'^{t) is obtained from the Hamiltonian ( |60| ) by 
replacing Xr{t) A,-(^)/(— 1)"^^*^. Its time-independent 
part Hq is given by (pQ) at t = 0, and the correction is 

ct)*. 



H'^{t) = A{cos(7r7V) - cos[Sxr{t) - ttN]}{c - 
{ck - c{)dk- 



* = 



The perturbation H'^{t) vanishes at T/Ec 
limit the spin contribution to the correlator (|5: 
the Green's function of operators U, 



(63) 
(64) 



0. In this 
becomes 



Xi°)(r) = (r,C/(r)C/(0))o, 



(65) 



where (. . .)o denotes averaging over the equilibrium ther- 
mal distribution with the Hamiltonian Hq. The explicit 
analytic result for this quantity is given by formula (66) 
of Ref. The result is an even function of the gate volt- 
age TV, and therefore within the approximation H'^{t) — 
the thermoelectric coefficient Gt vanishes. 

To find the leading contribution to Gt at small T/Ec 
we expand Eq. (|6^) to first order of the perturbation the- 
ory in H!^ (t) . The correction to Kg (r) has the form 



dt{TMt)U{T)U{0))o- 



(66) 



This corre ctio n is evaluated with logarithmic accuracy in 
Appendix |C 2|, 



X(i)(r) -^|r|2sin(27riV)ln 



Ec 
T + T 



T = 



SjEc 



\r\^cos^{7rN)- 



(67) 



(68) 
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It is important to note that although this result is the 
first-order correction in H'^{t), it is non-perturbative in 
the reflection amplitude r. 

Substituting the correlator K{t) in the form ( j^) with 
Kc and Kg given by Eqs. (|3|) and (]67)into the expres- 
sion for the thermoelectric coefficient (22|), we obtain 



DTT e he 1 + 1 

, [a;2-H(r/r)2]cosh2(x/2)' 



(69) 



At temperatures T ^ F this expression reproduces the 
perturbative result (|50|). The latter is valid until T ^T, 
and at T <C r the thermoelectric coefficient Gt becomes 



TT^ Gl T3 1 sin(7rA^) ^ 

(-rT — — — ~ 1 — ~ ~ TTT m ■ 



1 



6O72 



M I 1 2 cos3 (TTiV) I r 1 2 cos2 (tt A^) ■ 



(70) 



The dependence of Gt on the small reflection amplitude 
illustrates the non-perturbative nature of this result. It 
is also worth noting that at low temperatures the de- 
pendence of Gt on the gate voltage iV is strongly non- 
sinusoidal. 

To find the thermopowcr S — Gt/G one can use the 
expression (|69| ) and the non-perturbative result of Ref . |lj 
for the conductance G of the SET, 



G = 



GlT 
8jEc 



(x2 + n'^)dx 



[a;2 + (r/r)2] cosh^(a;/2) 



(71) 



At relatively high temperatures T ^ T the thermopower 
S — Gt/G obtained from Eqs. (|69|) and ( [7l| ) coincides 
with the perturbative expression ^5^) . In the more inter- 
esting case of low temperatures T <C F, we find 



TT^ 1 T 1 

5* — tan(7rA^) In — -r — -. 

5 eEc |r|2cos2(7rAr) 



(72) 



The new energy scale F arising from the non- 
perturbative solution is always small compared to the 
charging energy, see Eq. (|6|). It is important to keep in 
mind that F is a function of the gate voltage, and vanishes 
near the Coulomb blockade peaks, A^ = ±i,±|,±|,.... 
As a result, even at T <C Ec\r\'^ the perturbative re- 
sults (|^) and (|l|) are still valid near the conductance 
peaks, whereas in the valleys the new asymptotics ( |70| ) 
and ( |72| ) apply. The crossover between these asymptotics 
occurs at the values of N where F = T, i.e., according 
to Eq. (H) at a distance SN - ^T/^ckP from the 
centers of the conductance peaks. At these points the 
thermopower reaches its maximum absolute value S'max, 
which can be estimated by substituting N = ^ + SN in 
either Eq. ([5l|) or Eq. ([72|), resulting in 



Ec T 



(73) 



0.2 - 
to 0.1- 



T = 


0.3 E^. 


T- 


0.125 £p 


T = 


0.025 £p 


T- 


0.005 


T- 


0.001 E^ 




0.2 



0.4 0.6 

Gate voltage A' 



FIG. 2. Thermopower of SET as a function of gate voltage 
at different temperatures. The curves for 5* = Gt/G are ob- 
tained numerically from Eqs. (^) and (^l|) at \rf = 0.1 and 
T/Ec = 0.3, 0.125, 0.025, 0.005, 0.001. As the temperature 
is lowered, the amplitude of the thermopower oscillations first 
grows according to Eq. (^) and then decreases in agreement 
with Eq. (^3|). Below the crossover temperature Ec\r\^ the 
shape of the oscillations becomes non-sinusoidal. 



The exact shape of the Coulomb blockade oscillations 
of the thermopower found from Eqs. ( |69| ) and (71) is il- 
lustrated in Fig. 0. 



IV. SUMMARY AND DISCUSSION 



We presented a theory of the thermopower of single 
electron transistors in the regime when the coupling of 
the quantum dot to one of the leads is strong. The the- 
ory is applicable to devices with relatively large dots, 
where the effects of finite quantum level spacing can be 
neglected, and the main transport mechanism is inelastic 
cotunneling. Using the fact that the coupling to one of 
the two leads is weak, we obtained the expression ( [2^ ) 
for the thermoelectric coefficient Gt in terms of the cor- 
relator K{t) describing the charging of the dot strongly 
coupled to the other lead. The general expression ( [2^ ) 
is applicable to contacts with arbitrary coupling. We 
applied it to the case of coupling via a quantum point 
contact with a single transverse mode and almost perfect 
transmission, |r| <C 1. In the case of spin-polarized elec- 
trons we found sinusoidal Coulomb blockade oscillations 
of the thermopower with the amplitude ^ e~^\r\T/Ec, 
Eq. (^5|). Experimentally the polarization of electron 
spins can be achieved by applying a strong magnetic 
field. In the absence of the magnetic field the ther- 
mopower is given by the ratio of non-perturbative ex- 
pressions (p9) and (vw- At relatively high tempera- 



tures T ^ i?(7|r|2 the Coulomb blockade oscillations of S 
are sinusoidal, with the amplitude ~ e^^\r\'^ \n{Ec /T) , 
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Eq. (^T|). At lower temperatures T <^ Ec\r\^ the oscil- 
lations are non-sinusoidal, Fig. ^, and their amplitude is 
given by Eq. ([73|). 

We are aware of only one experiment on the ther- 
mopower of SET in the strong coupling regime, Ref. 
In this experiment the Coulomb blockade oscillations of 
the thermopower S{N) were measured at different val- 
ues of the reflection coefficient. Only one published 
curve S{N), measured at |rp = 0.2 ± 0.1 approached 
the strong tunneling limit \r\ <^ 1. In this case the ther- 
mopower remained sinusoidal even at the lowest avail- 
able temperatures. To observe the more interesting non- 
sinusoidal behavior of S{N) one would have to mea- 
sure the thermopower at lower temperature to reach the 
regime T ^ Ec\r\'^- This may require making a sam- 
ple with a larger quantum dot to ensure that the lowest 
temperature is still large compared to the quantum level 
spacing. 
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APPENDIX A: TUNNELING DENSITY OF 
STATES 



In this appendix we present the derivation of Eq. (|lj) 
for the tunneling density of states i^(e). We start with 
the standard expression for the density of states 



K6) 



2tt 



[G^(6)-G^(6)] 



(Al) 



where and G'^ are the retarded and advanced Green's 
functions, which can be obtained by the analytic contin- 
uation of the Matsubara Green's function G(e„). 

In the frequency representation the Matsubara Green's 
function can be written as 



G(e„) 



Jo 



dr exp(ie„T)G(r), 



(A2) 



FIG. 3. Representati on of the deformation of the 
T-integration contour in Eq. (A2). For positive e„ the contour 



should be distorted into the upper half plane, and for negative 
e„ — into the lower half plane. 

Keeping in mind that the retarded Green's function 
G^(je„) = G(e„) at e„ > 0, and using the fact that for 
the fermionic Matsubara frequencies exp(ie„/3) = —1, we 
can then express the retarded Green function as 



G (ie„) = i / dtcxp(-e„t) 
Jo 

x[G{it + 0) + G{f3 - + it)]. (A3) 

The analytic continuation to real frequencies can now be 
performed in the last line of Eq. ( [A^ ) through the sub- 
stitution e„ — > —ie. 

We then obtain a similar expression for the advanced 
Green's function G^{e) using the relation G'^(ze„) = 
G{en) at Cn < 0. Combining the two results, we find 
the following expression for the density of states (Al) 



1 

2^ 



dteicp{iet)Giit + 0) 

1 r°° 

— / dte^p{iet)G{P -O + it). (A4) 
2vr J_oo 



Since the Green's function G(t) is analytic everywhere 
except on the lines Rct = 0,±/3,±2/3, . . ., we can shift 



the integration contour in the first line of Eq. (A4) by 
t ^ t — if3/2 and in the second one hy t t + if3/2. As 
a result we obtain Eq. (p^. 



where e„ = 7rT(2n -\- 1) are the fermionic Matsubara fre- 
quencies. Depending on the sign of e„ the r integration 
contour can be distorted to the upper or lower half plane 
as shown in Fig. 



APPENDIX B: K{t) FOR SPINLESS ELECTRONS 

In this appendix we derive the results (|2|) and (^l]) for 
the correlator K{t) in the spinless case. 
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1. Evaluation of Ko{t), Eq. (28) 



To derive Eq. (^8|) we evaluate gaussian integral (|2£ 
under the assumption S' — 0. First, we find (l)r{x,t) 
that mini mizes the action So + iSc(t). Differentiating 
Eqs. (27a) and (|27b|) with respect to (j){x,t), we find 



d^(t)r + v%dl4>r - 2vfEc 



nr{t) H cf) 



N 



5{x) ^ 0. 



The solution of this equation has the form 



£;cexpf-^) 



(Bl) 



where ujn = 27rnT are bosonic Matsubara frequencies, 
and nT-{ujn) is the Fourier transform of ririt), Eq. (pj). 



n^(w„) 



1 



In the calculation of the correlator Ko{t) the integrals 
over the fluctuations of the field (p{x, t) about the saddle 
points (prix, t) and (jy^ix, t) in the numerator and the de- 
nominator of Eq. ( p5| ) cancel each other. Thus Ko{t) is 
given by the ratio of the saddle-point values of th e resp ec- 
tive integrals. Substituting (Bl) into (27a) and (27b) we 
find the saddle point action in the form 



[So+Sc{r)], 



Ec ^ 1 - cos(27rrTn) 



= 1 "("+2;F%) 



In the denominator of Eq. (|2^) we have the saddle 
point action at r = 0; according to Eq. ( p35| ) it vanishes. 
In the numerator of Eq. (25) the time r is finite. Assum- 
ing T > £'p^ and T < Ec, we find 



- [So+Sc{T)]4,=^^(^^t) ^ 21n 



^Ec \ sin(7rTT)| 



(B4) 



The corr elator Ko{t) is now found by exponentiation of 
Eq. (p|). The result is given by Eq. (||). 



2. Evaluation of K(t) to first order in r, Eq. (31) 



To derive the first-order correction (|T]) to the corre- 
lator Kq{t) one has to evaluate the gaussian functional 
integral (p3). It is convenient to integrate with respect 
to fluctuations 93 = of the field (j){x, t) about the 
saddle point (j)r{x,t). Then the integral ( pO[ ) takes the 
form 



{S')r = Re / 
Jo 



D. 



, (B5) 



where the averaging (. . .) is performed over the fluctu- 
ations around the saddle point <pr- This averaging can 
be viewed as integral (|3^) with Ut and N in the charg- 
ing action Sc set to zero. The evaluation of this inte- 
gral is straightforward, but lengthy. It can be avoided 
by noticing that the expression in the square brackets in 
Eq. ( B5 ) is time- independent and has the meaning of the 
first-order correction to the ground state energy of the 
Hamiltonian (23a)-(23c) at fi = iV = 0. Substituting its 
value found in Ref. |13|, we get 



{S')r 
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(B6) 



Using Eqs. (Bl) and (B2) we now find 



{S')r {S')o 



1 



r\Ec Re e 



2iTTN 



(B2) where 



X / di(e^[^(*)-^(*--)l-l), (B7) 



(B8) 



^ ' ^ ,T -X- -fee 



n=l " + 2^ 



At T ^ Ec the series can be evaluated explicitly for 
arbitrary i. 



Fit) 



_ f2g:cot(^Tt), 



2 Jo ^y- 



i{Ecty/Tr) 
1+y 



t:$>Ec\ 



(B9) 



In order to find the thermoelectric coefficient (^23) we 



need to find K(t) at r ~ > E^^. At these time 
scales the details of the short-time behavior of e*^'-*' are 
irrelevant, and one can replace 



1 



5{t) -l-isinF(t), 



where the constant f 

dx 



Ec 

« 1.59 is defined as 
siii{xy) 



1 - cos 2 



dy- 



l + y 



(BIO) 



(Bll) 



Substit utin g the approxim atio n ( |B10 ) into the integral 
in Eq. (|B7|) and using Eq. (p9|), we find 



(5%-(5')o-27e|r|[cos(2^iV) 



2tt^T 
Ec 



sin(27rA^) cot(7rTT)] . (B12) 



The calculation of the conductance G and the ther- 
moelectric coefficient Gt requires the knowledge of the 
even an d od d in r components of K{t), respectively. 
In Eq. (B12) we retained only the leading-order terms 
in T /Ec for each of these components. Substituting 
Eq. (pi|) into (||), we arrive at Eq. (IsT 
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APPENDIX C: K{t) FOR ELECTRONS WITH 
SPIN 



1. Evaluation of Kc(t,t';r) and K.s{t,t'), Eq. (42) 



The correlator Kc{t,t';T) defined by Eq. (12a) can be 
presented in the form 

Kc(i,i';r) = ^Re[Kt{t,t';T) + K-{t,t';T)], (CI) 



where 



(C2) 



The calculation of the correlators amounts to evalu- 
ation of gaussian integrals. Similarly to the calculations 
of Appendix B2, it is convenient to integrate over the 
fluctuations ipc about the saddle point 

0,^.(0, i) = ^ _ ^-^(^^n.(c.„)e— (C3) 



where n{ujn) is given by Eq. (B2h. The saddle point 
Eq. ( |C3D is easily obtained from ( p31| ) by replacing Ec ^ 
2Ec, nr{t) n^(t)/\/2, and N N/V^. Substituting 
(t)c{0,t) = (j)c,TiO,t) + (p{t) into Eq. (|c|) we find 

Kf {t, t'; t) = exp {zV2[0c,r (0, t) ± (/.,,,(0, t')]} 

X exp {~2{ipS)[Vc{t) ± iPcit'm . (C4) 



To evaluate the last factor in Eq. ( p4D we introduce the 
generating functional 



W[{J{uj,-,)}] - ( exp 



(C5) 



This gaussian integral is completely determined by the 
saddle point value if'l{t) of the field (pc, 



W[{J{LOn)}] = exp 



lT^J{uJn)ipi{-CJ„) 



Next we note that fluctuations of (pc{t) coincide with 
those of 0c (0, i) a,t N = rir — 0- Then Urit) in 
Eq. ( |38bD plays the role of a source term similar to 
J{t). More precisely, they are related according to 
J{t) = {2V2Ec/TT)nr{t). Then the saddle point ipi{t) 
can be determined from ( |C3| ) at iV = and ririuJu) = 
{TT/2y/2Ec)J{cJn), and we obtain 



W[{JM}] = exp 



J(w„)J(-a;„) 



2Ec 



(C6) 



Differentiating the functional W with respect to J{uon) 
and J{—ujm)i from Eqs. (03) and (|C^) we find 



(V3c(-Wri)</5c(Wm)) = 



2r 



2gc 



(C7) 



In the time representation this result takes the form 

TT giw„(t-t') 
(^c(%e(t')> = ^T^— — ^ 



(C8) 



The asymptotic behavior of this correlator is 

— r, \t~t'\«^, 

(C9) 



iln— , 



Substituting Eqs. (O), (p3), and (pq), into Eq. (^) we 
find the correlator K^t^ t^r) at |t — rp;» E^^ in the form 
(H), where 



in ( n + 



Be 



This definition of Xr(i) can be rewritten in the form (|4^). 

Our derivation of Kc(i, i'; t) allows one to find Ks{t, t') 
as well. Indeed, at Eq the actions of the charge 
and spin modes are identical. Taking the limit Ec — > 
in Eqs. (|C3|) and (^), from ^ and (^) we find the 
correlator Ks{t, t') at \t - t'\ > D^^ in the form (0). 



2. Evaluation of the correlator A'i^'(r), Eq. (|67D 



In this Appendix we outline the derivation of the cor- 
relator ( |67| ) starting from Eq. (|66|). At small temperature 
T <C Ec one can expand the expression ( |63| ) for i?^(t) 
to first order in (5%,- and present Eq. (pq) in the form 



$(r,i) = (rt(c + ct)t^'(i)C/(T)[/(0))o. 



(CIO) 
(Cll) 



Here we introduced the shorthand notation {c ~\- c^)t = 
c{t) + c\tY To evaluate $(r, <) we substitute the ex- 
pression (61) for U. Since the operator (c — c^) com- 
mutes with the Hamiltonian _ffo, the Green's function 
(Tt(c — c^)r(c — c^)o)o = —1, and we find 

$(t, t) = {Tt{c + c))t^{t){c + ct),(c + ct)o)o. (C12) 

Considering that the Hamiltonian is quadratic in 
fermion operators, one can use Wick's theorem and 
present ^(t, f) in terms of the single particle Green's 
functions: 

$(r, t) = Gi (r)G2 (0) + d (< - r)G2 (-i) 

-Gi(i)G2(T-t). (C13) 
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Here Gi{t) and G2{t) are defined as 

Gi{t) = {Tt{c + c^)t{c + c^)oh, 
G2{t) = {Tt{c + c^)t^{0))o. 



(C14a) 
(C14b) 



Evaluation of tlie Green's functions (C14) can be facili- 
tated by noticing that upon the substitution A cos{ttN) = 
A the Hamiltonian Hq coincides with the Hamiltonian 



H 



dk 



(C15) 



in Eq. (44) of.Ref. |T^. This Hamiltonian was diagonalized 
to the formtj 



H = E 



(k(ClCk + ClCk)dk, 



(C16) 



where E is the ground state energy of the Hamiltonian 
H and the fermion operators Gk and Ck are given by 



Ck - C^fc 



V2X 



-r2 
r 



V2 



d^k' Ck' - 



Ck = (cfc + cl 



(C17) 
(C18) 



Here F = AttX^/vf, which in our notations becomes 



To find t he Gr een's functions (C14) we invert the trans- 
formation (017) and obtain 



Ck 



V2 



r2 



dCk' d{a')Gk'+0hMcU, 

o{^k)Ck + e{-a)Gl, 



a- 



(C19) 



dk 



{Ck + Cl) . (C20) 



Using the se re sults, the definition of Eq. (64), and 
the form ( pl6 ) of the Hamiltonian, we easily obtain the 
Green's functions 



2r 

Gi(t) = — sgni 

TT 



VF J-, 



d^ 



^d^ 



r2 r 



^2 + r2 e/5« + 1 ■ 



(C21) 
(C22) 



The Green's functions Gi{t) and G2{t) are odd and 
even functions of t, resp ectively. Noticing also that 
Sxrit) given by Eq. ( 45b ) is invariant with respect to 
the change of variables t ^ t ~ t, we conclude th at the 
contributions of the second and third terms in Eq. ( |C13| ) 
to the integral (GIO) are equal to each other. Fi nally , 
the first term in Eq. ( C13 ) does not contribute to ( C10| ) 
because it is independent of i, and the ti me in tegral of 
Sxrit) vanishes. Therefore we rewrite Eq. (CIO) as 



iir^) (t-) ^ _2Asin(7riV) / dxTit)Gi{t)G2iT - t)dt. 
Jo 

(C23) 

Without loss of generality we can assume F ~ T ^ Ec- 
Since Sxrit) ^ -{■n'^T /2Ec)coi\TTTt] near t = 0,/3, 
see Eq. (11), the integral ( |C23| ) diverges logarithmically 
at t ^ and t ^ (3. These divergences are cut off 
at the short time scale E^^ and the long time scale 
min{F-i,T-i}. Due to the fact that Gi{+0)G2{t) = 
— G'i(/3— 0)6*2 (t—/3), the two divergences add up. There- 
fore with logarithmic accuracy the correlator ( G23 ) is 
given by 



27rA 



kP{t) = —sm{TrN)Gii+0)G2{T) In 



Er 



E, 



T + r 



(C24) 



Substituting G'i(-I-O) = 1 and the expression ( |C2^ ) for 
G2, we arrive at Eq. (67). 
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